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Abstract 



We perform a semiclassical calculation of the magnetoresistance of spin- 
less two-dimensional fermions in a long-range correlated random magnetic 
field. In the regime relevant for the problem of the half filled Landau level 
the perturbative Born approximation fails and we develop a new method of 
solving the Boltzmann equation beyond the relaxation time approximation. 
In absence of interactions, electron density modulations, in-plane fields, and 
Fermi surface anisotropy we obtain a quadratic negative magnetoresistance 
in the weak field limit. 
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The problem of two-dimensional transport in spatially random static magnetic fields 
(RMF) attracted a lot of attention over last few years. There exists a number of experimen- 
tal realizations of two-dimensional (2D) RMF created by randomly pinned flux vortices in 
a type-II superconducting gate, grains of a type-I superconductor, or a demagnetized per- 
manent magnet placed in the nearby of the 2D electron gas |I| . A formally similar problem 
arises in the contexts of the gauge theory of high Tc cuprates and the composite fermion 
theory of the half filled Landau level (HFLL) ||, §. 

Various analytical and numerical results obtained so far seem to indicate that in the 
case of a continuum system with unbound spectrum all states are localized |^ while in the 
lattice case there might be an extended state at the center of the band 0. The effective low- 
energy description in terms of the unitary cr-model constructed in [0] suggests that the RMF 
problem belongs to the unitary random ensemble. Compared to the conventional problem of 
potential scattering corresponding to the orthogonal case, the effects of quantum interference 
in the RMF are suppressed as a result of broken time-reversal symmetry. In particular, the 
logarithmic temperature dependent weak localization corrections to the conductivity appear 
only in the next order in the metallicity parameter kpl » 1: 6wi(J = ~ ^2fc^; log ^ (where 
the temperature dependence comes from the inelastic phase breaking time T^p{T) in the 
regime t^{T) » Ttr = ;^). It suggests a larger localization length Lioc ~ I exp{^{kpiy) 
than in the orthogonal case. 

It was shown in IQ that quantum fluctuations of particle's positions in the direction 
transverse to their classical trajectories are strongly reduced by the RMF. And so is the 
contribution of quantum backscattering which leads to a negative magnetoresistance (MR) 
in presence of a weak uniform external magnetic field B. 

Therefore, in the RMF problem one might expect that the semiclassical transport the- 
ory is applicable in a wider range of length scales and temperatures than in the case of 
random time-reversal invariant potential scattering. Yet, a semiclassical treatment of the 
RMF scattering beyond the relaxation time approximation may lead to a non-trivial MR. 
This classical contribution due to the bending of particle's trajectories dominates over the 
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suppressed effects of weak localization, at least, at not very low temperatures. 

The semiclassical approach to the RMF problem was undertaken in a number of publi- 
cations. In the Subnikov-de Haas (SdH) oscillations of pxx{B) in the strong field limit 
were studied by summing over classical cyclotron orbits. However, to find 
Pxx{B) in the weak field limit before the onset of the SdH oscillations one has to develop a 
semiclassical analysis in terms of particle-hole pairs rather than single fermions. 

The authors of [|l^ used the linear Boltzmann equation where the RMF played a role 
of a random driving force instead of including it into the collision integral. The condition 
which makes it possible is kp^ » 1, where ^ is a correlation length (scale of a typical 
spatial variation) of the RMF. However the analysis performed in [0] was restricted onto 
the case relevant for the experiments [|l|] where, first, strong potential scattering provides a 
large bare value of p^x, and, second, the spatial correlation of the RMF has a finite range 

In the present letter we address the case when pxx is solely due to the scattering by the 
RMF described by the (not necessarily short-range) white- noise correlator < bgb^q >= f{q) 
independent of the external field B. We will not restrict our consideration onto the lowest 
order of perturbation theory. Therefore we will be able to discuss the case /(g) ~ e"^''^ 
relevant for the problem of HFLL where the spatial correlation decays as and the lowest 
order result simply diverges. 

We start with the Bolzmann equation for the distribution function g{t, r, 0, e) = go{^) + 
6g{t, r, (p, e) which reads as 

(«^ + z{;^V + «(fie + ^)^)<5(7(t,r,0,e) =2et;^i^, (0.1) 

where = vpn^ is a vector of Fermi velocity normal to the circular Fermi surface param- 
eterized by the angular variable 0, E stands for an infinitesimal external electric field, and 
qq = 6{6f — e) is the unperturbed Fermi-Dirac distribution. By contrast to the case of ordi- 
nary potential scattering (which can be also treated at fc^^ >> 1 as a random electric field) 
the Eq.(l) implies a trivial dependence of 6g on e, namely, Sg(t, r, 0, e) = —Sgit, r, (f))6{ep—e). 



We note, in passing, that one can use the Eq.(l) even in presence of long-range retarded 
gauge interactions between fermions which spoil the Fermi liquid coherence. Even though 
there are no well-defined fermionic quasiparticles, the kinetic equation in terms of local 
displacements of the fiuctuating Fermi surface can still be derived |1T[, thereby making a 



close contact with the idea of bosonization of 2D non-Fermi liquids |]T^ . 



The perturbed distribution function can be found in terms of the (retarded) Green func- 
tion of the Bolzmann operator for a given RMF configuration 

Sg{t, r, ^, e) = iej dcf)' J df^G{t; r, r"; 4>')v^,E5{eF - e) (0.2) 

and then averaged over all possible configurations. As usual, the averaging over disorder 
restores the translational invariance in both r- and 0-spaces: << G(t; r, r^; 0, 0') >>= 
Q(t; r — r'; (j) — 0'). Then the calculation of the frequency-dependent conductivity amounts 
to computing the p-wave harmonics of the g* = component of the Fourier transform of 

I 27r 
where i?c = ff" is the cyclotron radius of fermions with density Ue'- 
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In absence of the RMF the bare Green function is given by its harmonics gI^\u!) = — , 
In our case of no potential disorder {6 = 0^) the lowest order (bosonic) self-energy 



correction to Qi ^(cu) = cu — IVlc + T.i(uj) found in |T0| makes it only possible to study the 



5-dependent correction to the bare p^x at high frequencies or magnetic fields. To proceed 
with a more complete account of the effects of the RMF we first solve the equation for the 
Green function 

(z|- + iv^V + z(fie + ^-^)^))G{t - t'- r, r'; 0, 0') = 5{t - t')6{f- r')6{4> - 0') (0.3) 
at m 0(p 

for an arbitrary RMF configuration by Fourier transforming with respect to t — t' and f—r'. 
The solution of the resulting equation 



{to - v^q + iv^V + iiSl^ + ^)—))G{uj- q, f; 0, (/)') = 5(0 - (/)') (0.4) 

m acp 

can be searched in the form 

G(f, q- <p, <p') = -I dr j (^0"e^^«o '(«1</>.'^")e^*{r;r=-;0",<^') (q.S) 

Here ^ can be viewed as a "eikonal" phase of a particle-hole pair propagating in the RMF. 
Substituting (5) into (4) and integrating by parts we obtain the equation for ^: 

(,|. + V)* = -e-*e-^-^o ^M^l^e'-^o^e^* (0.6) 
OT m d(p 

Provided the RMF correlation length is large enough {kp^ » 1) one can linearize the 
Eq.(6) (we will comment on this point below) and then end up with an explicit solution 

*(t; f; 0, 00 = 1 £ dr'b(r- 4')^<^(0 - <l>'), (0-7) 

-* ^ _d_ _^ 

where — RcC " ^'I'n^ is a classical trajectory corresponding to the Larmor precession 
along the cyclotron orbit. Now taking the gaussian average of the exponent e** in (7) over 
different realizations of the REM we obtain the averaged Green function 

g{r,Q-<p-4>') = -i rdTe^^('^+^"=*)exp(i f dn f dT,F{2R, sin V ""'^ )^)5{<t^ ' <f>') 
Jo 2 Jo Jo 2 0(p^ 

(0.8) 

where F{r) is a Fourier transform of the RMF correlator f{q). The use of the Eq. (8) yields 
the conductivity 

(^xxi'^) — -I dr COS (Q,c + <^)t dr'(r — r')F(2RcSm—^)) + (u ^ —u) 

m Jo Jo 2 

(0.9) 

The formula for axyico) differs from (9) by the additional factor tein{flc^u!)T in the integrand. 

The exponential factor in the integrand in (9) is controlled by the amplitude of the RMF 
correlator F{r). To check the validity of (9) in the perturbativc regime we expand the 
exponent up to the first order in F{r) and perform the r integration first (it has to be done 
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with an infinitesimal exponent e ^'^ inserted into the integrand which specifies the retarded 
nature of the Green function (8)). Then one can readily obtain 

CTxxi,^) = l/rtri^c ~Uj)~'^ + {UJ ^ -Uj) 

where the RMF scattering rate is given by the expression 

roo J- 

l/Ttr(B)= dr COS QcT F (2Rc sin (0.10) 
Jo 2 

which is nothing but the result of the first Born approximation (BA) ||lO|. One can also 
obtain the formula (10) by estimating the ~ r term in the exponent in (9) in the r — oo 
limit which would be equivalent to the relaxation time approximation. 

In what follows we will concentrate onto the case relevant for the problem of HFLL 
corresponding to f{q) = An'^ariee''^'^^ where we introduced the dimensionless coupling a 
proportional to the density rii of ionized Coulomb impurities separated by the spacer of the 
width ^ from the 2D electron gas and the amount of gauge fiux quanta $ attached to every 
electron [a = ^^"^^ for a filling factor u with an even denominator $ Q]). In the picture 
of composite fermions (CF) 0, ^ each charged impurity also becomes a source of a gauge 
magnetic fiux. It is believed that scattering by randomly distributed fiuxes provides the main 
mechanism of the CF momentum relaxation while the potential scattering is negligible. 

At -B = the Eq. (10) gives the CF elastic transport rate = ||a which coincides 

with the BA result found in [Q. Although this estimate certainly holds for a << 1, it is no 
longer valid at large a when the CF mean free path (MFP) Icf = 2^a^^ gets shorter than ^. 
The case of HFLL (z/ = 1/2) appears to be marginal {a = 2) and / determined this way just 
equals ^. Since a is proportional to $^ the situation becomes even worse for compressible 
states at fractions with higher even denominators. 

This observation signals about a failure of the BA at large coupling a when the matrix 
element Mp-_p» ~ describing a single event of CF scattering by a typical magnetic 

impurity satisfies neither of the two conditions M^^^ << or required for the 

validity of the BA [|T^. This is the semiclassical regime of a small (A0 ~ << 1) angle 
scattering, which, however, can not be treated in the lowest order of the perturbation theory. 
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Nevertheless, if the parameter kp^ is large enough, so that 

kp^ » v^, (0.11) 

one can resort on the so-called eikonal approach [|T^ which was essentially implemented by 
the above solution of the Boltzmann equation. In fact, it is the condition (11) which allows 
one to solve the equation (6) in the first order in 6(r^. 
The use of (9) then leads to 

2 2 

= ^{kpl) = £-(fc^e)e"i^i(a), (0.12) 

where Ki{x) is the modified Bessel function of the second kind. At small a the Eq.(12) 
reproduces the above BA result while at strong coupling it predicts the MFP to be / ^ ^ \/^- 

Thus the condition (11) actually provides that kpl » 1 which is necessary to verify the 
very use of the kinetic equation (1) and to obtain the metallic value of the conductivity. 

In the intermediate coupling regime relevant for the case of HFLL the condition (11) is 
fairly well met {a = 2 and kp^ ^ 15) and the Eq.(12) gives the CF conductivity a^f = 
^(kFC)e'^Ki{2) which is 2.06 times greater than the BA result. It is worthwhile mentioning 
that the experimentally measured resistivity at z/ = 1/2 is about 3 times smaller than the 
BA estimate 

To add to this point we mention that a somewhat similar problem arises in a semiclassical 
calculation of the amplitude of SdH oscilattions of the MR at large deviations from v = 1/2 
0. The calculation carried out in reproduces the experimentally observed scaling ^£2^2^. ^ 
exp{— const/ B^jj^) with the effective field -Be// = B — Aime » m/r^^^ but gives a factor in 
the exponent which is about 2300 times greater than the measured one. Presumably, such 
a discrepancy also stems from the inadequacy of the lowest order approximation (which 
systematically overestimates the strength of the RMF scattering) and can be fixed by going 
beyond it. 

One could wonder if the strong coupling behavior of the MFP as a function of the strength 
of RMF fluctuations (/ ~ ^ rather than / ~ ^) can be already seen in the self-consistent 
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BA which in the case of a finite-range potential scattering is apphcable at ,^ << 5"^/^. 
Adapting this method to the RMF problem one can achieve a self-consistent improvement 
of the lowest order result by inserting a factor e^'^^'^"' into the integrand in (10) and solving 
the resulting nonlinear equation for Ttr- 

For the case of HFLL this equation reads as 

Vp JO Vp 

and gives the MFP l{a) in agreement with the above eikonal calculation. In terms of a 
conventional diagrammatics the equation (13) corresponds to the sum of uncrossed "rain- 
bow" diagrams for the bosonic self-energy Si. By contrast, the eikonal result (8) includes 
contributions of some of the crossed diagrams in all orders of the perturbation theory. 

The MR in the RMF problem was identified in as the -B-dependence of Ttr- It is 
obvious, however, that in absence of other mechanisms of momentum relaxation one can not 
use the Eq.(lO) to analyze the behavior of the MR at Vtf. < l/r^^. 

Moreover, in the HFLL case the long-range character of F{r) ~ r^^ leads to a logarithmic 
divergency of the second derivative of Pxx{B) at i? — 0, which is, of course, an artifact of 
the above expansion. In the framework of the more accurate self-consistent BA one obtains 
that in the weak field limit Ttr{B) is an increasing function of B which implies a negative 
MR. Our self-consistent analysis shows that the decreasing behavior of Ttr{B) (positive MR) 
in the case of /(g) ~ e"^'?^?^ and small a reported in [10|, in fact, holds only for large enough 
fields: Qc > l/^tri^^TT^- 

The more complete eikonal treatment of the RMF problem confirms this prediction. 
Expanding (9) up to the second order in Q^. we obtain the MR defined as ^^""^^"^^ = 

pxx 



XX 



+ 1 < r dr'(3r'V - 4r'^)F(r') >) (0.14) 
< 1 > 12 JO 



p^^ 2 < 1 > 

where 

<t" >= drr"exp(- / dr'ir - t')F{t')) 
Jo Jo 



Notice that in the relaxation time approximation (< >= /q°° (irr'^e"'^/'^*'') the expression 
(14) contains only two terms proportional to < > and < r >^ which exactly cancel out 
and result in zero MR. Since the exponential factor Jq dr'^r — t')F{t') appearing in all our 
calculations instead of r/rtr behaves as ~ at small r, the combined effect of these terms 
on the MR is negative and can dominate over the remaining (strictly positive) contribution. 

The MR remains negative for all couplings in a wide class of realistic RMF correlation 
functions /(g) including the HFLL case at a = 2 when the Eq.(14) yields 



2 

rj 



Pxx 

while in the strong coupling limit of large a the MR becomes ^p^'^^^') = (^]^rtr)^(f — !)• 

The possibility to apply our conclusions to real systems is comphcated by such factors 
as Fermi surface and fermion dispersion anisotropy, periodic electron density modulations, 
or even a finite tilt of an external field which are all known to yield positive contributions 
to the experimentally measured MR. 

The existing artificial realizations of the RMF all have a property that the random field 
depends on the applied external one. A theoretical analysis shows that the RMF correlation 
function f{q) appropriate for the references (a) and (c) in |I| is proportional to B while 
the one relevant for the reference (b) is ~ B^. In either case this trivial -B-dependence 
necessarily leads to a positive MR and does not allow one to explore a more sophisticated 
origin of a non-zero MR which is the only one present in the case of f{q) independent of B. 
The quadratic positive MR was also reported in a numerical simulation of the lattice 



version of the RMF modeled as an assembly of uncorrelated random fluxes |jT^. Again, the 
anisotropy of the lattice fermion dispersion and the obvious lack of the large parameter kp^ 
in the short-range case preclude us from making a direct comparison with our results. 

In the case of HFLL one may also think of the observed positive MR as resulting from the 
interference between the RMF scattering produced by charged impurities and strong gauge 



interactions of CF. It was recently shown that this interference might be responsible 
for the strong non-universal logT correction to Pxx{T) at u = 1/2 |jl6|]. The account of the 



gauge interactions drastically changes the renormalization properties of the RMF problem 
and can even reduce the localization length down to Lioc ~ /exp(j^||^). 

The zero field analogue of such an interference phenomenon (the Altshuler-Aronov cor- 
rection) is known to lead in some cases to a positive MR [0. In our case of interest the 
effect of the external magnetic field could be even more pronounced. The work intended to 
clarify this issue is in progress. 

To summarize, we develop a new method of solving the semiclassical Boltzmann equation 
for the RMF problem beyond the relaxation time approximation provided kpl » 1. The 
obtained solution is used to calculate the RMF transport time, mean free path, and zero- 
field conductivity which all appear to be greater than the corresponding results of the Born 
approximation. In particular, we propose a new estimate of the semiclassical conductivity 
of composite fermions in the compressible state at u = 1/2 which is about twice the value 
found in [Q] and agrees better with the experimental data. We also show that in the ideal 
case of non-interacting fermions with uniform density and isotropic circular Fermi surface 
the magnetoresistance in the RMF is negative and quadratic in the weak external field. 

The author is indebted to Per Hedegard for a valuable discussion of the results of this 
paper. 
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